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Abstract 

Stroboscopic wavepacket basis sets [P. Bokes, F. Corsetti, R. W. Godby, Phys. Rev. Lett. 101, 
046402 (2008)] are specifically tailored for a description of time-dependent processes in extended 
systems like non-periodic geometries of various contacts consisting of solids and molecules. The 
explanation of the construction of such a basis for two simple finite systems is followed by a 
review of the general theory for extended systems with continuous spectrum. The latter is further 
elaborated with the introduction of the interaction representation which takes the full advantage 
of the time-dynamics built into the basis. The formalism is applied to a semi-analytical example 
of electronic transport through resonant tunnelling barrier in ID. Through the time-dependent 
generalisation of the Landauer formula given in terms of the Fourier expansion of the transmission 
amplitude we analyze the temporal character of the onset of the steady-state. Various time-scales 
in this process are shown to be directly related to the energetic structure of the resonant barrier. 
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I. INTRODUCTION 



Time-dependent density functional theory (TDDFT) allows the formulation of first- 
principles simulations of many-electron problems in the language of non-interacting par- 
ticles moving in the effective time- dependent Kohn-Sham potential. Actual solution of such 
a problem requires a particular representation, or a basis set, for one-electron wavef unctions. 
There are many possibilities to choose from: system independent basis sets like the real- 
space representation [X] and plane-waves 2, or more-or-less system-specific basis sets like 
various types of localized orbitals, either of Gaussian type J4] or numerical orbitals 5), Q], 
ground state canonical Kohc.Sha.n orbrtals correspond to the studred system, normally 
localized Wanier functions [7| and possibly many other. The system independent basis set 
have the advantage of being general and not preferential to any particular system, which is 
balanced by the problem that typically one needs very large number of basis functions to 
obtain converged results. This has been the main motivation in the development of basis 
sets that are less general, more specialized for the description of the involved physics or 
chemistry at sufficient precision. With the latter the almost ubiquitous criterion was their 
description of the ground states properties of the studied system. 

In this paper we discuss a new type of basis sets that are specifically tuned for a description 
of time-dependent processes in extended systems like solids or non-periodic geometries of 
various contacts consisting of solids and molecules f|. The major complications for these 
systems is their size. In the ground state one may hope to use "nearsightedness" of the 
electrons and study the relevant questions using cluster or periodic models. In non- 
stationary processes this is not possible in general since electrons may traverse long distances 
and 'see' arbitrarily remote parts of the system. A problem related to this, and one that 
is particularly relevant for the description of quantum transport within the TDDFT is the 
possibility that the exact, or at least sufficiently precise exchange-correlation potential might 
need to be very nonlocal [k]] even though the numerical importance of this effect in not clear 



at present 
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One possible way is to study very large systems and look only at times smaller than the 
time it takes for the electrons to realise that their box is in fact finite , lU] . Another 
possibility how to treat infinite systems is to employ embedding self-energies that 
play the role of semi-infinite parts of the total system, often called electrodes, assumed 
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to behave in a known way, and the rest, typically a small finite subsystem that can be 
treated fully within the TDDFT calculation. The stroboscopic wavepacket basis set gives 
a framework that avoids the use of embedding self-energies in principle and still treat the 
system as infinite. Individual basis functions are tailored to the time- dependence in the 
electrodes and the time-propagation of each electron there is analytical. Only in the regions 
where the system differs from the electrodes it is necessary to propagate the occupied orbitals 
numerically. Still, even here the wavepacket basis gives a very appealing picture of electrons 
being redistributed in between different wavepackets due to local disturbances in the effective 
potential. 

In the following section we introduce the concept of stroboscopic construction of a basis 
set on a very simple examples: two level system consisting of Is and 2s orbitals of Li 
atom and a subset of the eigenstates of Harmonic oscillator. In Sec. Ill Bl we review the 
formalism of the wavepacket basis for any Hamiltonian with a continuum spectrum which 
is in Sec. Ill CI further elaborated into interaction representation based on the wavepacket 
basis. The following two sections give an example of the use of the basis on transport of 
electrons through a simple resonant potential barrier. We find a simple and direct relation 
between the Fourier analysis of a transmission amplitude and the temporal behaviour of the 
relaxation into the current-carrying steady state. 



II. STROBOSCOPIC WAVEPACKET BASIS 



A. Finite systems 

Let us consider a two level system, e.g. the canonical Is and 2s orbitals obtained within 
the ground state DFT calculation of the Lithium atom, | Is) and |2s) with eigenenergies e\ 
and e 2 . To obtain the ground state density, we occupy the first orbital twice and the second 
one once with the result 

n(r) = 2| (r\ls) \ 2 + | (f\2s) | 2 . (2.1) 

In principle, any unitary rotation between these occupied canonical orbitals gives identical 
density. We will keep one electron, say with its spin down, in the Is orbital, but the other 
two, both with spin up, will be placed into new orbitals, 

\g) = ^(|ls) + |2s» (2.2) 
3 



\u') = ±={\U)-\28)) (2.3) 
This is, of course, a well known fact. New and interesting point of view arises if we view the 



state \u') as the result of time-propagation of the |o)[29l| 



e 



|^ = J_ e -ieir + e -i(e 2 - ei )r ^ = ^ie.r |^ for T = * _ ^ 

a/2 e2 — ei 



This means, we could consider orbitals \g) and |w) = e~ tHT \g) as two orthogonal states which 
can be filled with electrons. This rather trivial example clearly demonstrates the principal 
idea of the stroboscopic basis set: to obtain orthogonal basis by means of a propagation 
from suitably chosen initial state. 

The second example where such a construction works for discreet spectra is harmonic 
oscillator. Consider N consecutive states with eigenenergies e n = ftw n, n = n , . . . ,n + 
N — 1 and eigenstates \4> n )- We first choose the initial state as the equal combination of the 
eigenstates, 

1 no+N-l 



n=no 

In principle, the phases of each of the coefficient in the initial state are arbitrary, which 
reflects certain freedom in choosing a particular basis. 

The remaining members of the new basis are obtained through propagating the initial 
state with time step r = 2tt/(Nljq), 

\i> m ) = e~ lilmT |^o) , m = 0,l,...,iV-l. (2.6) 

It is a simple excercise to show that this set is orthogonal, (V'mlV'm) = &mm'i an d spans 
identical subspace as the chosen eigenstates. The new basis carries certain resemblance with 



the coherent states used extensively within the quantum optics [17j. Namely, both lead 
to almost classical interpretation of the dynamics of a single quantum particle, e.g. the 
contribution of a single electron to the total density is oscillating from right to the left of 
the parabolic potential with the period T = 2tt/ujq. On the other hand, in contrast with the 
coherent states, the time-propagated states by construction consists only from eigenstates 
with energies within given interval which allows for exact description of non-interacting 
many-electron ground states or simple non-equilibrium models, as it is described in the 
following sections. 
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B. Infinite systems 



Systems of infinite extent in space have not only discreet but most importantly also 
continuous spectrum of eigenvalues. In our previous work [8] we have shown that the states 
with eigenvalues from chosen interval of energies (e, e + Ae) in the continuous spectrum can 
be unitarily transformed into wavepacket basis, in a similar way to the examples in the 
previous section. In the following we will review the basic steps of this formalism. 

Let the system of noninteracting particles under consideration be initially characterised 
by a reference Hamiltonian H with a continuous part of its spectrum covered with intervals 
of energies (e n ,e n+ i), n = 0,1,... to which we refer here as the energy bands. We will 
assume that each energy within given band is N n times degenerate. The eigenstates of the 
reference Hamiltonian, 

H |e, a) = e |e, a) , 
need to be taken normalized to the delta-function of their energies, 

(e',a > \e,a} = 5(e-e')8 a , a ,. (2.7) 

From the above set of eigenstates we can generate the wave-packet basis set (WPB) by 
first choosing the initial set of wave-packets (WPs) 



£n+l 



|n,0,a) = — [=/ de'\e',a), ra = 0,1,2,... (2.8) 



for each energy band {(e n , e n+ i)} n and each value of the degeneracy index a. Ae n = e n+1 — e n 
represents the width of all the considered energy bands. 

The construction of the WPB is completed by forward and backward time propagation 
of each member of the initial set 

\n,m,a) = e- iAmTn \n,0,a) , m = ±l,±2,... (2.9) 

by regular, band-dependent time steps r n = 2n / Ae n . This choice of time steps guarantees 
orthonormality of consecutive wave-packets within each band since, 

(n, m, a\n, m', a) = -r— f dede' (e,a\e~ t{e ' m '~ em)T \e' ,a) 



A 



e. 



I dee- l < m '-^ r = 5 m>m ., (2.10) 
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where we have used the delta-function normalisation (I2.7p . The same equation also guar- 
antees that WPs with different degeneracy index a, as well as WPs from disjunct energy 
bands are mutually orthogonal. Collecting together all such WPs together with the discreet 
eigenstates (if they are present) of the reference Hamiltonian we form a complete orthogonal 
set of states. The regular time steps between individual WPs remind a view through a 
stroboscope on a continually moving single wavepacket. This analogy led us to refer to such 
a complete orthogonal set as the stroboscopic wavepacket basis. 

A particular example of the stroboscopic basis are the orthogonal wavepackets construc- 
tured from the plane waves and used first by Th. Martin and R. Landauer 18| . In fact, 
this work served as a motivation for generalisation into the stroboscopic construction in the 
present form. 

There are several very attractive properties of this basis that we would like to emphasize: 
(1) instead of continuous and delta-function normalised states we can work with countable 
even though infinitely many states that are localized in space and normalised to one; (2) 
occupying all states from all energy bands below the Fermi energy, Ep, we recover identical 
many-particle ground state of non-interacting particles characterised by the reference Hamil- 
tonian; (3) while the total density of such many-particle ground state is time-independent as 
it should, electron in each WP state in the energy interval (e n , e n + Ae n ) can be interpreted 
as making transition from m-th WP to m + 1-th WP in time r n . As it will be demonstrated 
in Section IHIl all these properties can be very useful for quick and insightful description of 
non-equilibrium processes in many-electron systems. 



C. Time- dependent wavepacket basis 

In the previous section the wavepacket basis set was constructed in terms of stroboscopic 
snapshots of time-propagation of initial set of suitably chosen wavepackets. Clearly, this 
construction is valid for any initial time. This can be exploited even further within the in- 
teraction representation where each time-dependent basis function is an evolving wavepacket, 
and different orthogonal basis functions are time-shifted with respect to each other by r. 
This means that any wavefunction <j)(x, t) = (x\<f), t) occupied by an electron can be expanded 
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as 

(x\(f), t) =} j (x\i, a, n; t) c ij(X>n {t) (2.11) 

i,n,a 

where i runs over different energy bands, a through different degenerate states within the 
band and n over different orthogonal wavepackets from the given band i. Due to their very 
definition through time-evolution these time-dependent basis function have the property 

(x\i, a,n;t + Tj) = (x\i, a, n + 1; t) , (2-12) 

for any real time t, T{ given by the width of the energy band i and any integer n indexing 
the various WPs within the band. 

There is not much gained from using this particular representation if the actual Hamil- 
tonian of the system is the reference Hamiltonian which generates the wavepacket basis at 
all times. In such a case the coefficients Cj )Q , jn (t) are either equal 1, if the considered band 
i is below the Fermi energy, or equal 0, if the band i is above the Fermi energy. However, 
when the actual Hamiltonian locally differs from the reference Hamiltonian at times t > to, 
i.e. the Hamiltonian is 

H = H + V(t) (2.13) 

where V{t) ^ only in finite part of the overall system, Q C -R 3 , and t > to, the wavefunc- 
tions' coefficients will evolve according the Schroedinger equation in the interaction picture 

d 

ifaCj,a,m(t) = 22 (j,a,m;t\V(t)\i,a',n;t)c ij0l/tn (t). (2.14) 

i,a',n 

Since the perturbation V(t) is assumed to be nonzero only in Q, and the wavepackets are 
localised in space, the above matrix elements are nonzero only for those WPs that have 
significant amplitude in Q. The coefficient Ci^ n (t) for times before the corresponding basis 
function (x\i, a, n; t) is evolved by the reference Hamiltonian into the region Q will be equal 
1, reflecting the ground state for t < to imposed by the reference Hamiltonian in the past. 
During the following time, when the matrix elements of this WP basis function are nonzero, 
the coefficient c i)Q , jn (t) will change according to Eq. (I2.14p . Finally, once the WP basis 
function (x\i, a, n; t) departs from Q, the coefficient will freeze and will be kept constant 
once again. 
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It is clear from this discussion that the use of the interaction picture representation 
of the stroboscopic wavepacket basis allows us to translate the time-dependent quantum- 
mechanical problem of a system with infinite extent into one that demands numerical solution 
of finite number of dynamical equations (12.141) . We will demonstrate this strategy in the 
following sections where we address the behaviour of the relaxation into a current carrying 
steady-state for a model resonant nano-junction at moderate-to-long time-scales. 



III. ELECTRONIC TRANSPORT THROUGH RESONANT TUNNELLING BAR- 
RIER 



A. Definition of the model 



For the purposes of demonstration of the use of the time-dependent stroboscopic basis we 
will study a very simple model of resonant tunnelling. In many respects, such a model can be 
viewed as a prototype of many interesting phenomena observed in the nanoscopic transport 
through molecules attached to metallic electrodes 
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21] . The time-dependence in 



tunnelling 



e.g. 
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las been addressed using various wavepackets by many authors in the past, 
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271 ] . The distinguishing contribution of the present treatment is 



the correct many-electron occupation corresponding to finite bias and low temperatures and 
hence its direct consequences for the current-voltage characteristics. The system will be 
strictly one dimensional, with two delta-function potential barriers localised at x = and 
x = a with equal strengths given by a parameter A (see Fig. [I]) . 

Applying voltage to such a model results in decreasing the constant potential to the 
right of the double barrier by a value —V. The value of the potential in between the two 
potential barriers is chosen to be —V/2. The physics of this choice is that the electrons 
inside the double barrier will effectively screen the field and the total drop in potential will 
be symmetrically split into the two barriers where the electronic density is low and hence 
the screening is weaker. 

The simplicity of this model allows for analytical construction of its right- and left- going 



S 



scattering eigenstates for any positive energy e, e.g. for the right-going state we have 



e ikx _|_ r (^ & ikx^ _ X < 

k \x) = < C 1 (e)e^ + C 2 (e)e-^, q = ^2(e + V/2) < x < a , (3.15) 

t(e)e iKX 1 K= y/2(e + V) x>a 

where all the functions r(e), Ci(e), C 2 (e) and t(e) can be obtained from the conditions on 
continuity the wavefunction and finite discontinuity of its derivative at x = and a, and are 
simple functions of the energy e of the considered state. Only the transmission amplitude, 
t(e), is directly relevant for the present work and its functional form is 

^ (k + q + 2i\)(k + q + 2i\)e i{ - K -^ a - (k - q + 2i\)(k - q + 2iX)e^ K+ ^ a ' 



The prefactor of the scattering state (p§(x), 1/V2nk, guarantees correct normalisation given 
by Eq. (12.71) . The right- and left- going scattering states at the same energy are examples 
of two degenerate continuum eigenstates. In the language of the general formulation of the 
wavepacket basis [Eqs. (12.81) . (I2.9P ], their indices R and L represent two possibilities of the 
degeneracy index a. The resulting WPs for a = R and L will be arriving to the barrier from 
the left or right respectively. 

The particular values of the parameters will be chosen in the following way: the Fermi 
energy is set to Ep = 1.0, the strength of the barriers will be varied through values A = 
0.1, 1.0 and 2.0. We will use two different distances between the barriers, a\ = 9.5 and 
Gt2 = 10.5 to clearly demonstrate the different types of relaxation behaviour. Both of these 



are above the Fermi wavelength X F = 2n / (y/2Ep) = AAA to allow for existence of several 
resonances within the energy interval (0, Ep) at zero bias. In Fig. [2] we show the transmission 
of such a potential for nonzero bias V and several values of A. 

B. Non-linear response to abrupt switching-on of the voltage 

In equilibrium, all from the right- and left- moving wavepackets arising from the energy 
bands below the Fermi energy will be occupied so that overall current is zero. This is not 
the case in non-equilibrium situation, and nonzero current can be obtained. The calculation 
of the current alone can be done using the free electron wavepackets (FWP), 

" en+1 de 



{x\i,a,n;t)= I -^ = e iakx e~ iet , k = \[2~t (3.17) 



i.e. the wavepackets where the reference Hamiltonian is that of free electrons Hq = 
— (1/2) d 2 /dx 2 . This Hamiltonian characterises the regions to the right and left of the dou- 
ble barrier. (We will use the round bracket to distinguish the FWP from a general WP 
\i,a,n;t).) The degeneracy index a for FWPs is equal to + or — and indicates the right- 
and left- moving FWPs respectively. Within this representation the current due to the 
energy band or the bias window (E F — V, E F ) (the corresponding time step between neigh- 
bouring FWPs is Ty = 2n /V) is just a simple counting: number of electrons that occupy 
n-th right-going FWP divided by the time it takes for an electron to move from n-th to the 
n + 1-th FWP, 

= NM = m v> (3,8) 

Ty Z7T 

where N n (t) is the average occupation of the n-th FWP. 

Immediately after the switching-on the bias, the evaluation of the current at a chosen 
position of the FWP and time needs to be analyzed numerically since the time-dependent 
Hamiltonian will occupy various unoccupied bands of higher energy and these would also 
contribute to the current. However, it is possible to arrive at a simple expression for the 
occupations N n (t) for moderate to long time-scales when only WPs in the bias window 
(Ep — V, Ef) contribute to the current. Under these conditions the average occupation is 
given in terms of the overlaps between the scattering wavepackets correspondingly the final 
Hamiltonian which includes the barrier and the bias potential, and the FWPs [8], 

— oo 

N n (t) =2^\(n,+,i \i,R,m;t) | 2 , (3.19) 

m=0 

where 



1 



(i, +,n \i, R,m;t) = — / det(e) exp{— ie(t + (m — n)ry)} (3.20) 
V Je f -v 

The z-th energy band is understood to be the bias window, the factor 2 stands for spin 
degeneracy. The Eq. f 1 3 . 2 j) is strictly valid only for intermediate and long times times, 
t > 1/Ep, as well as only for FWPs that are sufficiently distant from the resonant barriers, 
n > a/iypTy). The simplicity of Eqs. (13.191) and fl3.2Qj) is striking. The Eq. (I3.20j) is simply 
a Fourier analysis of the transmission on the interval of energies that fall within the bias 
window. Let us introduce the coefficients of this Fourier series as 

tv(l) = J yt(e)e-^exp{-ie(lry)}, (3.21) 
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where the time is split into the integer multiple of the time step ry and the rest, t = 
N t Tv + < i] < Ty. Using this convention we can write the total current at the position 
of n-th FWP as 

— oo 

I n (t) = ^\t r ,(Nt + m-n)\ 2 , (3.22) 

m=0 

which can be readily evaluated. It turns out that this sum converges rapidly: for the here 
studied resonant barrier the results for the times show in the Fig. H] and [5] are converged if 
we take just 15-20 terms. This is an indication that the number of needed WPs to achieve 
convergence of results does not need to be large. Similar observation has been found in 
convergence of the ground states density for square potential barrier [28( . 

We now proceed to discuss the results for two different values of applied voltage. First 
we choose such a voltage that only one resonance appears in the bias window (Ep, Ep — V). 
Next we increase the voltage and shorten the distance between the two barriers to = 9.5 
to capture two resonances. In the long time limit, the current goes into its steady-state 
value Jqo given by the Landauer formula 

-Ep 

de\t(e)\ 2 . (3.23) 

E F -V 



The resulting I — V characteristic for our model is given in the Fig. [3j Whenever new 
resonance enters the bias window (Ep — V, Ep) a rapid increase in the current appears, 



which is a well known feature in the transport through molecular nano-j unctions 20]. 

The resulting dependence of the current on time, studied at the location of the n = 4-th 
FWP (n = corresponds to FWP located behind the position x = a) for the bias voltage 
V = 0.3 is shown in Fig. HJ In this case only a single resonance enters the bias window, as it 
is indicated in the inset of the Fig. |H It takes time t = 4ry for the front with nonzero current 
to appear at the 4-th FWP basis function in the right lead, and the steady-state current is set 
in with damped oscillations whose period is given by the time step Ty. This is particularly 
clear for transmissive enough barriers (A = 0.1). Such a behaviour is not just an artefact 
of the WP basis and the choice of the energy band of width V; identical voltage-dependent 
oscillations were found using non-equilibrium Green's function formalism [15]. Increasing the 
strength of the barriers, i.e. increasing the parameter lambda, the oscillations are suppressed 
and instead, new time-scale appears in the functional form of the time-dependence: a slow 
exponential approach towards the steady-state current. The time-constant of this approach 
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is given by the width of the resonance present within the bias window V. This can be 
confirmed by inspection of the transmission spectrum (inset of the Fig. H]). 

If two resonances are present within the bias window (V = 0.7, Fig. ED, the behaviour 
of the current relaxation is similar to the previous case for highly transmissive barriers. 
However, increasing A and hence making the barriers more opaque brings in a new time- 
scale: the slow exponential approach to the steady-state value determined by the widths of 
the resonances is now superposed with oscillations with period slightly larger than Ty. This 
can be ascribed to the energy difference between the two resonances that are both within 
the bias window (the inset of the Fig. [5]) . 

The example of electronic transport through resonant barrier shows how the time- 
dependent expression for the current in terms of the Fourier decomposition of the transmis- 
sion, Eq. (13. 2 1 p allows for a very simple analysis of the time-dependent processes involved 
in the relaxation of the current towards its steady-state value. 

IV. CONCLUSIONS 

The stroboscopic wavepacket basis represent a novel basis set specifically tailored for a 
description of time- dependent processes in extended systems We have explained the idea of 
the stroboscopic construction on a simple two level system and a finite system consisting 
of subset of eigenstates of the harmonic oscillator. In the following we have reviewed the 
general theory for construction of the stroboscopic wavepacket basis for an extended system 
with continuous spectrum. We have introduced interaction representation within this basis 
which practically reduces the number of dynamical equations that need to be solved to finite 
number. 

We have applied the time-dependent formalism to semi-analytical example of electronic 
transport through resonant tunnelling barrier in ID. The use of stroboscopic basis allows 
for physically appealing and at the same time mathematically exact formulation of the 
transport. Furthermore, through the time-dependent generalisation of the Landauer formula 
in terms of Fourier analysis of the transmission amplitude we could identify the temporal 
behaviour of the establishment of the steady state with nonzero current where various time- 
scales could be put directly into relation with the energetic structure of the resonant barrier. 
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Fig. 1 Form of the model resonant potential at nonzero bias V. The value of the po- 
tential inside the barrier is chosen to be halfway between the bias. 

Fig. 2 Transmission function for the model resonant potential at zero bias. Ep = 1.0 
guarantees that several resonances are occupied in equilibrium, before the bias is 
switched on. 

Fig. 3 The steady-state current-voltage characteristics of the resonant barrier. Rapid 
increases in current correspond to an emergence of next resonance within the 
bias window. The green arrows indicate voltages for which we have studied the 
time-dependent approach to the here indicated values of current. 

Fig. 4 Relaxation of the current towards its steady-states value, 1^, for the case when 
single resonance is present within the bias window. Varying the strength of the 
barriers alters the character for oscillatory to exponential. The inset shows the 
transmission spectrum at this value of the bias with the resonance present in the 
bias window. 

Fig. 5 Relaxation of the current towards its steady-states value, 1^, for the case when 
two resonances are present within the bias window. Varying the strength of the 
barriers alters the character for oscillatory to exponential with superimposed os- 
cillations that can be ascribed to the energy distance between the two resonances. 
The inset shows the transmission spectrum at this value of the bias and the two 
resonances present in the bias window. 
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FIG. 1: Form of the model resonant potential at nonzero bias V . The value of the potential inside 
the barrier is chosen to be halfway between the bias. 
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FIG. 2: Transmission function for the model resonant potential at zero bias. Ef = 1.0 guarantees 
that several resonances are occupied in equilibrium, before the bias is switched on. 
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FIG. 3: The steady-state current-voltage characteristics of the resonant barrier. Rapid increases in 
current correspond to an emergence of next resonance within the bias window. The green arrows 
indicate voltages for which we have studied the time-dependent approach to the here indicated 
values of current. 
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FIG. 4: Relaxation of the current towards its steady-states value, 1^, for the case when single res- 
onance is present within the bias window. Varying the strength of the barriers alters the character 
for oscillatory to exponential. The inset shows the transmission spectrum at this value of the bias 
with the resonance present in the bias window. 
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FIG. 5: Relaxation of the current towards its steady-states value, 1^, for the case when two 
resonances are present within the bias window. Varying the strength of the barriers alters the 
character for oscillatory to exponential with superimposed oscillations that can be ascribed to the 
energy distance between the two resonances. The inset shows the transmission spectrum at this 
value of the bias and the two resonances present in the bias window. 



19 



